Derealization transition in an exactly solvable many-body system in Two Dimensions 
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We study derealization transition in a many body system in two dimension. We identify the 
presence of a complex vector potential that gives rise to this transition. 
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Non-Hermitian Hamiltonians appear in the study of 
quantum dots coupled to leads where it is convenient 
to exclude the degrees of freedom in the leads by in- 
tegrating over them. This scheme results in an effective 
non-Hcrmitian Hamiltonian whose eigenstates contain an 
imaginary part that accounts for the decay into the leads. 

In Ref. [2| , Hatano and Nelson, considered a different 
kind of non-Hermitian Hamiltonians, one with a constant 
imaginary vector potential. This leads to the appear- 
ance of a first order derivative in the Hamiltonian that, 
as pointed out by Efetov |lj , can be assimilated to the 
appearance of a certain preferential direction in space. 
The introduction of this ingredient gives rise to new ef- 
fects such as a localization transition in a disordered one 
dimensional system 

In this context, we study analytically an exactly solv- 
able iV-body system in two-dimensions, which shows a lo- 
calization to derealization transition even in the absence 
of random potential. The hamiltonian is non-hermitian 
with a complex vector potential. 

We derive the ground state eigenvalue and eigenfunc- 
tion for such systems and extract several interesting 
properties by identifying the square of the radial wave 
function with the joint probability distribution (JPD) of 
eigenvalues of the Laguerre ensembles of random matrix 
theory (RMT). We also study the particle density and 
pair correlation function (PCF) for such systems. We 
observe "universality" in the PCF as expected for such 
systems from random matrix results. 

Let us consider a two-dimensional system of N par- 
ticles in the x — y plane, subject to a vector potential 
A = [xx + yy]/[2(x 2 +y 2 )}- The single body Hamiltonian 
is given by 



H sp — 



{Ti/i)V + iA 



The many body Hamiltonian is given by 



H = Y,m/i)v+ l A ] ], 



(i) 



(2) 
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where Aj is the vector potential. 

Making the transformation x = r cos 6 and y = r sin ( 
and taking K = m = 1, we get 



Here, we have used the relation 

d 2 Id 2 Id 



V 2 = 



dr 2 r 2 dO 2 r dr 



(3) 



(4) 



In terms of the polar coordinate, the vector potential can 
be written as 



A = -jv, 

f 

2r~ 



where, 



and 



ld_ 

rdO 1 



xcosO + ysmU, 



V r = 



dr ' 



(5) 



(6) 



= —x cos + y sin 0, (7) 



where x, y, f and 9 are unit vectors along the x, y, r and 9 
directions respectively. This expression in turn allows us 
to view the complex vector potential as some velocity op- 
erator acting along the angular coordinate. This maybe 
the reason for which one can do away with the random 
part of the potential (as is seen in the work of previous 
authors), but retaining the effect of disorder. However, 
one needs to understand this aspect more clearly. 
Replacing this in the Schrodinger equation 



HV n (r i ,6 i )=E n y n (r i ,6 i ), 



(8) 



we obtain the ground state eigenfunction and eigenvalue. 

Let us take the solution wavefunction of the form 
y n {ri,0i) = R(n, . . . ,r N )@(6i, . . .,9 N ). Using the sep- 
aration of variable and equating the angular part to zero, 
we get 
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^ T^R{ri, . . . ,rjv) = E n R(n, ...,r N ) (9) 



The angular part of the equation is given by 



(10) 



Now, let us assume that the solution wavefunction has 
the form 



Rn(n) = (r< - rj) A exp(-a fc i 



±i ■ 



e(^,...,M=^exp(— ^0 



(11) 



(12) 



where A can take only two values, i.e. and 1. Then 
A = will correspond to bosons while A = 1 to fcrmions. 
In that case, ak = J2i K ~F ik n ] and the eigenvalue for the 
TV-body Hamiltonian will be given by 



E, 



-N(k 2 



(13) 



We must note that for a given system k is constant. The 
ground state (the state that minimizes the real part of 
E n ) is localized for k n = 0. The localized ground state 
is Eq = —Nk 2 . In the extended state, in order to have a 
well defined wave function, we need to specify the bound- 
ary conditions specific to the system. The quantity K is 
a measure of the amount of localization. 

Now, we will prove that (|llf> is indeed a solution of the 
radial equation ©. Let us take 



IK' 



ip = exp(-a fe ^2 '• 



(14) 
(15) 



part of the Hamiltonian, A can only be or 1. Thus we 
have proved Eq.^T]. The proof of the angular part is 
straightforward. 

It is interesting to note from Eci. (|13|) that for k n = 0, 
the wavefunction is exponentially localized and the eigen- 
value does not contain any imaginary part. However, as 
k n increases, the wavefunction becomes extended in na- 
ture and a complex part is added to the eigenvalue. 

It is at this point that, we will write ip n as 

i>n = C 1/2 J] \n - r 3 \ X exp[±a k n] exp[±^], (20) 



Defining the scalar product {ipni^Pm) as 

(Vn,^m) = C J J Y[\n - r/ X 

cxp(-K y^(fj) ± i(k n - k m ) ^ ' 

— 7T CS nm . 



(21) 



Here for n ^ m, the term on the right is highly oscillatory 
and hence goes to zero. For n — m, we get the normaliza- 
tion. Here, C is the well known Selsberg constant given 
by 



N-l 



c- 1 = n r(1 + A +; ?A) ? +jA) . (22) 



3=0 



r(i + A) 



Then \ip n \ 2 is given by 



\Vn\ 



\n~r/ exp(-2 K ^(r l ), (23) 



i<3 



with P = 2A, and C being the normalization constant. 

Now one may interpret \"4> n \ 2 to be identical with the 
j.p.d. of Laguerre ensembles of random matrices Q. We 
define the m-particle correlation function 



It can be easily seen that 



, d 2 ip d(p dip 



(h 



dn dr.; 



Substituting Eqs.(14,15), the first term gives 

,0V 



dr 2 



while 



i<j 



:<j)(p. 



(16) 
(17) 

(18) 
(19) 



The middle term becomes zero due to symmetry. We 
can also see that for Eq. ^T] to be a solution of the radial 



Y m (ri,r 2 , 



dr 



m+l 



. drivdO 



.dd N \tp \ 
(24) 



as the probability of finding m-particles in the intervals 
ri and rj -I- Ar^ and 9i + A9i , irrespective of the position 
of the other particles, m = 1 and 2 correspond to the 
particle density and PCF respectively. 

It has also been shown in Ref. [j] that Yi(x) corre- 
sponds to the density of zeros of the polynomial having 
weight function given by that of the Laguerre Polyno- 
mial. It has been shown by Dyson in the context of ran- 
dom matrices that for f3 = 1, 2 and 4, Y m can be written 
in terms of orthogonal and skew-orthogonal polynomials. 
For = 2, corresponding to the Fermionic solution, the 
PCF can be written as 



Y 2 (n,r 2 ) 



N-l 
fj,=0 



^MriM^Jexp^n] (25) 



3 



where (r) are orthogonal polynomials corresponding to 
the normalization condition 



qMqAv) exp[-2nr]dy = h^5^ v . (26) 



It should be noted that our choice of circular coordi- 
nate has considerably simplified our calculation as |V>| 2 
does not contain 9 explicitly. Using the Christoffel- 
Darboux formula in Ea.|25). we can perform the sum. 
Finally using the asymptotic form of the Laguerre poly- 
nomials, and taking the limit, we get for the particle 
density. Here, we will give an alternate derivation (used 
previously in Q ) for the derivation of particle density. 

We define the resolvent 



(27) 



which satisfies 



G{x + iO) 



Yi(r) 



dr- wrYi(r). (28) 



Also, from the definition of \ip n \ j we may write Q 

dYl{r) -/?/^r 2 -2 K y 1 (r). (29) 



Or 



(r - r 2 ) 



For large N, replacing Y2(r, r%) — Y\(r)Yx{r2) and drop- 
ping dYi(r)/dr, we get 



(z-r) 2 



(z - r)(z - r 2 ) 



(30) 



This gives 

f3zG 2 {z) + A K N -4kZG(z) = 0, (31) 
solving which we get for the particle density 



7T/3Yi(r) = 2k 



Np/n-r 



(32) 



It is interesting to note that as pointed out earlier, k is 
a measure of the degree of localization. For k large com- 
pared to N, particles get more strongly localized around 
the origin. However, as k decreases, particles start get- 
ting delocalized, until they form extended states. Also, 
it should be noted that in the thermodynamic limit, the 
particle density, after proper scaling, can be made inde- 
pendent of (3 and hence applicable for both boson and 
fermion. However, in this case, the nature of localization 
is exponential for bosons (A = 0), while it is given by 
Eq.[S21 for fermions (A = 1). 

To calculate the PCF, essential for the study of ther- 
modynamic properties, for /3 = 2, it has been proved in 
Ref.jij that for interparticle spacing {r\ — r 2 ) = Ar — > 0, 
N — > oo and defining r' — Ar.Yi(r), the scaled PCF is 
universal. For (3 — 2, it is given by 



Y 2 (r 1 ,r 2 ) = 1- 



sin 2 (7rr') 



(33) 



(wr') 2 ■ 

This shows that the probability of finding two particles 
in the same state (Ar = 0) is zero. This is in agreement 
with the Pauli exclusion principal, applicable for a system 
of fermions. 



Thus we have studied the localization to derealization 
transition in a system of A^-particles in the presence of 
a complex vector potential. For localized state, particles 
tend to gather around the origin, as shown in Ea.|32j. 
At the same time, the PCF shows that the particles are 
not allowed to condense to one position. The correspond- 
ing eigenvalue is real. However, for delocalized state, we 
observe that the eigenvalue has complex part. 



S.G. is grateful to Dr. Louis Foe Torres for many useful 
discussions. 
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